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Abstract 



The structural properties of a linear polymer and its evolution in time have a strong 
bearing on its anisotropic stress response. The mean-square bond length and mean bond 
angle are the critical parameters that influence the time-varying stress developed in the 
polymer. The bond length distribution along the chain is uniform without any abrupt 
changes at the ends. Among the externally set parameters such as density, temperature, 
strain rate, and chain length, the density as well as the chain length of the polymer 
have a significant effect on the stress. At high density values, changes in mean-square 
bond length dominates over changes in radius of gyration and end-to-end length. In other 
words, bond deformations dominate as opposed to changes in size and shape. Also, there 
is a large change in the mean-square bond length that is reflected as a jump in the stress. 
Beyond a particular value of the chain length, n = 50, called the entanglement length, 
stress-response is found to have distinctly different behavior which we attribute to the 
entanglement effects. Short chain polymers more or less behave like rigid molecules. There 
is no significant change in their internal structure when loaded. Further, temperature and 
rate of loading have a very mild effect on the stress. Besides these new results, we can now 
explain well known polymeric mechanical behavior under dynamic loading from the point 
of view of the evolution of the molecular dynamics and the derived structural properties. 
This could possibly lead to polymer synthesis with desired mechanical behavior. 
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Nomenclature 



Roman Symbols 

hi Bond vector between the united atoms i — 1 and i 

fij force on united atom 

G H*H 

H matrix consisting of edge vectors of the simulation cell 

I identity tensor 

Pij momentum of the united atom j in chain i in real coordinates 

Tij position vector of atom j in chain i 

Sij scaled coordinate of atom j in chain i 

f number of degrees of freedom 

Qk eigcn values of the mass distribution tensor, A; = 1, 2, 3 

Gxy element of mass distribution tensor 

K stiffness parameter in FENE potential 

ks Boltzmann's constant 

K^i bond torsion stiffness parameters, / = 0, 1, 2, 3 

K$ bond bending stiffness 

Li cell dimension in the direction e;, / = 1, 2, 3. 

niij mass of the atom j in chain i, In this case rn,j = m 

N number of chains 

n number of atoms in a chain 

R end to end length of the polymer chain 

Rq maximum allowable bond length 

Tc cut-off distance for short range potentials 

Rg radius of gyration of the polymer chain 

s dynamical variable to control temperature 

T temperature 

t real time 

Text external temperature imposed on the system 

U potential energy 

V volume of the simulation cell 

Wg inertia parameter corresponding to the extended variable s 
Greek Symbols 

Uij scaled momenta of atom j in chain i 

T virial stress tensor 

e energy parameter in Lennard-Jones potential 

A stretch ratio 

(j) dihedral angle 

ip chain angle, angle between end to end vector of chain and loading direction 

a length parameter in Lennard-Jones potential 

T scaled time 

6 bond angle 

^0 mean bond angle 

Subscripts 

i index of polymer chain, i = 1, 2, . . . A'^ 

j index of atom in a polymer chain, j = 1, 2, . . . n 

I coordinate directions, Z = 1, 2, 3 

Other Symbols 

( )' derivative with respect to scaled time r 

( ) derivative with respect to real time t 

^ Lagrangian 
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1. Introduction 

Polymers exhibit a wide variety of properties and behavior under different environ- 
mental conditions. Polymers are very large molecules made-up of smaller repeating units 
called monomers. Due to their flexibility in assuming random configurations they are 
capable of undergoing nonlinear and large elastic deformations when subjected to ex- 
ternal loads. Also, because of the relaxation in internal structure they show viscoelastic 
behavior as well. Further, these materials can be doped with foreign particles to increase 
their strength and stiffness as well as exhibit field sensitive properties. 

At present many models exist to describe the behavior of polymers that can be broadly 
classified as continuum models and molecular models. In continuum models there are 
those that are based on primitive mechanical elements such as mechanical springs and 
dampers, and are used in different combination to describe the behavior of polymers 
Wineman & Rajagopal (2000). Other continuum models are based on the free energy of 
the system, where an arbitrary form of the free energy is chosen to fit the experimental 
data, and the behavior of the polymer is described in terms of this functional form of the 
free energy using the laws of thermodynamics Thien (2002); Kankanala & Triantafyllidis 
(2003); Brigaduov & Dorfmann (2003). The molecular models of polymers are too many. 
A few that come to mind are the freely jointed chain model, phantom model. Rouse model, 
Zimm model, reptation model, network model, tube model, work-like chain model, etc 
Doi (k Adwards (I98G). None of these models describe the complete behavior of polymers 
in general but are only applicable in some limiting cases because of the simplifying 
assumptions that go into the model framework. Further, many of these models suffer from 
assumptions that could be physically unjustifiable. For instance, the phantom network 
model does not have any constraints on bonds crossing each other and is not physically 
plausible. Many of the advanced statistical mechanics models have been used to describe 
the behavior of polymers James (1947); Boggs (1952); Freed (I97f). Studies have been 
performed using some of these models and their ranges of validity established Gao & 
Weiner (1992); Rickayzen & Powles (2003); Kremcr et al. (1988). 

Molecular dynamic (MD) and Monte Carlo simulations too have helped in understand- 
ing the deformation behavior of polymers. Gao and Weiner Gao ^ Weiner (1984, 1987, 
19896, 1994) used a freely jointed monomer chain model of the polymer and came to the 
conclusion that stress anisotropy is primarily an effect of excluded volume interactions. 
The covalent bond force in freely-jointed models of a polymer network contributes pre- 
dominantly to pressure. As the volume fraction increases, contribution of covalent bond 
interactions to the stress anisotropy decreases Gao & Weiner (1989a, 1991a, &); Weiner 
& Gao (1994). 

Bower & Weiner (2004) studied the role of pressure on the mechanical behavior of 
polymers. They show using molecular simulation on an idealized elastomer, a relation 
between the difference stress and pressure. Difference stress is that which resists extension 
or stretch and pressure resists volume change. If the scaled density is less than one the 
difference stress is purely entropic. When the normalized density exceeds this value, 
then at lower temperatures difference stress contains an additional energy component. 
The role of monomer packing fraction on rubber elasticity was studied by Bower & 
Weiner (2006). They demonstrate that non-bonded interactions play a central role in 
controlling the mechanical response of rubber-like solids even under standard conditions. 
In their simulation, glass transition does not occur for scaled density less than 1, while 
at a scaled density value 1.2, they observe glass transition. They find that stress levels 
were considerably higher for higher densities as compared to lower densities. Finally, 
they conclude that short range excluded volume interaction is an active deformation 
mechanism contributing to difference stress whereas long range attractive interactions 
contribute only to mean stress or pressure. 

Weiner & Gao (1989) introduced the concept of intrinsic chain stress that generalizes 
the concept of axial entropic force for the case of inter- and intra-chain non-covalent 
interactions for dense systems. They show that the macroscopic stress tensor in a network 
can be expressed as a sum of individual chain stress contributions. They also find that the 
interchain non-covalent contribution to stress is non-hydrostatic in the deformed network 
and varies with deformation due to variation in chain orientation. They again use the 
concept of intrinsic stress to study the Mooney effect in swollen networks where they show 
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that due to non-hydrostatic contribution to stress, because of interchain non-covalent 
interaction, model polymer exhibits strain softening or Mooney effect characteristic of 
dry rubber solids Weiner & Gao (1990). 

Molecular dynamics simulation of polymer entanglement was performed by Krcmcr & 
Grcst (1990). They observed that the Rouse model provides an excellent description of 
short chain polymers while the dynamics of long chains can be described by the reptation 
model. They map the polymer to the experimental polymeric liquids and determine the 
entanglement length uniquely. They analyze the motion of the primitive chain to show 
the possibility of visualizing the confinement of motion of the chain within the tube. 

Stress anisotropy due to excluded volume interaction is quite well established in litera- 
ture on polymeric stress-strain behavior, the previous works as mentioned earlier. In the 
present work, we address the effect of internal structure on stress anisotropy in a polymer 
melt under dynamic strain. A generic polymeric system is modeled as an NVT ensemble. 
We study the stress-strain response of the polymeric system in combination with var- 
ious micro-structural parameters, namely, mean-square bond length, mean bond angle, 
mean-square end-to-end length, mean-square radius of gyration, mass-ratios, mean chain 
angle. The mean-square bond length and mean bond angle have significant effect on the 
stress anisotropy. Besides, bond deformation and conformational change becomes more 
significant when the chains are strongly aligned in the loading direction. With reference 
to external parameters such as density, temperature, chain length and rate of loading, 
we observe that the stress levels are highly sensitive to density and chain length. In con- 
trast, the temperature and loading rate are only mildly influential. These influences are 
correlated in terms of the micro-structure changes in the polymeric system. 



2. Mathematical model 

We consider a system consisting of N polymer chains each having n united atoms. The 
extended variables of the system are {s^-, s}, where, Sij is the scaled position vector of the 
atom j in chain i with respect to simulation cell dimension, s is an extended dynamical 
variable introduced in the Lagrangian to control the temperature Nose (1983). 



where 



H 



Hs, 



ii 
2.2 
L3 



(2.1) 



Alternately 



= Lis^ji, I = 1,2,3 (no sum on I). 



(2.2) 



Here Z = 1, 2, 3 refers to coordinate directions ei, 62, 63 respectively. Vij is the position 
of united atom j in chain i. Thus, Siji, Sij2 and Sijs S [0, 1), and the dimensions of the 
cell in the directions ei, 62 and 63 are Li, L2 and L3, respectively. The dynamics of the 
system is controlled by the Lagrangian 



N 



U{r,,) + ^--fkBT\ns, (2.3) 
i=i j=i ^ ^ 

which, when written in terms of scaled coordinates and scaled time, reads 



N n 

i=i j=i 



(s^.H*Hsy-C/(Hs,,) 



- fkBTlns. 



(2.4) 



In the above relations, / denotes the number of degree of freedom of the system. ( ) and 
( )' refers to the derivative with respect to the real time t and scaled time r {dr — dt/s) 
respectively. The momentum of the particle in terms of the Lagrangian is given by 
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(2.5) 



where G = H*H and 11,%,- denotes the scaled momenta. 

The equation of motion of a particle in a system whose Lagrangian is is given by 

dt dcii dqi 

Therefore, the equations of motion of the particles in the scaled coordinates can be 
written as 



and 



// 1-1 r -111 

^ii = — 2^ fjj ~ 2— Sj - G G Sj , 

J S 



(2.7a) 



WsS 



N 



1=1 j=l 



The equations of motion in real variables and real time now take the form 

771 S 771 



(2.75) 



(2.8a) 



s 

Ws 



N n 2 

-i=i j=i 



(2.85) 



The momentum, p^- , of the atom j in chain i in term of real coordinates and real time 
is given by 



H-*n, 



= m{vij — HH 



(2.9) 



The expression for the particle momentum in this ensemble is not the same as that 
in the NVE ensemble. There are some extra terms in the above expression because of 
variation in the cell dimensions. In the case when the system is held at a constant strain, 
its Hamiltonian remains conserved throughout the simulation. But when the system is 
subject to a time-varying strain, this is no longer valid. This accounts for the extra 
terms. While applying the strain on the system, all atom positions are scaled in the 
same proportion as that of the cell dimensions. Here we assume the material to be 
incompressible. So, if stretch ratio in the x direction is A, then the stretch ratio in y and 
z direction would be ^ because of symmetry. 

All the physical properties are defined in terms of momentum of the atoms thus defined. 
For example the kinetic energy is now defined as p^/2m instead of ^mr^. 

The calculation of stress is done by using the virial theorem Grciiicr (■/: al. (1997) of 
statistical mechanics. For a system of atoms interacting with two body potentials, the 
expression for the stress is given by 



1 
V 



iNkBT)I-(Y, 



dU 



\a>l3 



\r\a0 d\r\ai3 



Yaji ® ra/3 



a,/3 = 1,2,--- ,7Vn, (2.10) 



where r is the virial stress tensor, /c^ is the Boltzmann constant, T is the temperature, I 
is unit tensor and V is the current volume of the simulation cell. Various other parameters 
representing internal degrees of freedom and overall shape of the polymer chain considered 
in this study are briefly defined below. 

The mean-square end-to-end length of the polymer chains give information on the 
compact or open configuration acquired by chains on an average. It is defined as 
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(i?2) = (|r„-rip), (2.11) 

where r„ and ri are the position vectors of the last and first monomer or united atom, 
respectively, in a chain. 

The radius of gyration of a polymer chain denotes the entire mass distribution of the 
chain and plays a central role in interpreting light scattering and viscosity measurements. 
Radius of gyration for a polymer chain of length n is given by 



(^^) = ^(E|r^-f|') (2-12) 

In any configuration of the chain, the mean spatial distribution of the chain mass need 
not be spherical. The moment of inertia of the monomer mass distribution in a sense 
signifies the shape of the chain. Elements of the tensor describing the mass distribution 
have the form 



Gxy — — '^y^^l'ix ~ Tx){'''iy ~ fy)- (2-13) 

Eigenvectors of the above tensor give the principle directions in which mass is dis- 
tributed and the corresponding eigenvalues are a measure of the distribution. If the 
eigenvalues are arranged in a decreasing order, the two mass ratios of interest are the 
ratios of the second and third eigenvalues with respect to the first eigenvalue. These 
ratios give the relative measure of the mass distribution in the two transverse directions 
to the strain loading direction. 

Chain angle is the angle which an end-to-end vector of the chain makes with the 
direction of the loading. This is also a measure of alignment of the chains in the loading 
direction. It is given by 



where, e/ is the loading direction and / = 1, 2, 3. 

In MD simulations, the equations of motion. Equation 2.8a, need to be integrated to 
determine the new positions of the united atoms. The integration method that is quite 
popular in use is the Verlet algorithm Ra.pa.port (2(J()4). However, in this simulation since 
we are constraining the temperature, the accelerations are no longer only functions of 
position but also functions of momenta. This imposes considerable difficulty in integrating 
the equations of motion. A variant of the Verlet algorithm known as the velocity Verlet 
algorithm is used in this case since we need the current particle velocities too. This 
algorithm is very similar to the Verlet algorithm except for the fact that velocity rather 
than position at the previous time step is used to integrate the equations of motion. Due 
to application of temperature control, equations of motion of the particles depend on 
the current value of s and its first derivative. Again, the second derivative of s depends 
on the current particle momenta. This makes the equations of motion interdependent. 
This prevents the direct application of the velocity Verlet algorithm for integrating the 
equations of motion. Fox & Andersen (1984) gave a procedure to modify the velocity 
Verlet algorithm to overcome this difficulty without affecting the accuracy of the method. 
In this procedure, s and its derivatives are approximated by using the position and 
momenta at the previous time step. These approximate values are then used to integrate 
the equation of motion of the particles. The new values of the position and momenta 
thus obtained are used to make better approximation of s and its derivative which are 
used in the next time step. We also use the same method to integrate the equations of 
motion. 

We outline the principle numerical integration steps in this procedure. 



x{t + dt) = x{t) + dixit) + -{dt)'^x{t), 



(2.15a) 
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(2.15&) 



X \ t + ^dt) = i{t) + ^dtx{t), 



x(t + dt) = X 



t+ -dt 



-dtx{t 



dt). 



(2.15c) 



Now the equations of motion are integrated in the foUowing manner. 

(1) Assume that the values r{t),r{t), r{t) and s{t) are known at time t. With these 
values, evaluate r{t + dt) and s{t + dt) using Equation (2.15a). 

(2) The velocities r{t + ^dt) and s{t + ^dt) are next determined in the next half time- 
step using Equation (2.156). 

(3) Using Equation (2.86), determine s{t+dt) using the s{t+Y), s{t+dt) and momenta 
Pij{t). Since we are using the value of the temperature control variable s, its derivative 
and linear momenta from the previous time-steps to obtain s(t + dt), we refer to 
this value as Sapp{t + dt). This value in substituted in Equation (2.15c) to get the 
approximate value Sapp{t + dt). 

(4) Since the position of the particles is known at step t + dt, the force on each particle 
can be calculated from the gradient of the potential energy at this instant. 

(5) From Equation (2.8a) we observe that acceleration at time t + dt is a. function of 
current particle velocity. Further, from Equation (2.15c), particle velocity at time 
t -I- is dependent on the acceleration at the current time step. This implies that 
acceleration and momenta are interdependent. Substituting the expression of accel- 
erations from Equation (2.8a) into Equation (2.15c) we get 



Tij {t + dt) 



dt.s 



dt 
2" 



HH 



dt 



HH 



HH^ (HH^)* 



HH^ (HH-i)* ) HH~i 



dt 



dt iij 

2 TO 



(2.16) 



Substituting this expression for the velocity at the current time in Equation (2.8a), 
the acceleration at the current time step is obtained as 



Tij {t + dt) = 



TO 



dt s 

'2 s 



HH 



dt 



HH 



(HH 



HH 



HH 



I+fHH-i-(HH-i)*) If,, (t 



(HH 

dt 
+ 2 



HH 



(2.17) 



When the deformation is symmetric, HH ^ — (HH ^)* vanishes. In this case, the 
above expressions change to 



rij(t + dt) 



dt s 



dt 
*+2 



dt iij 

T TO 



dt 



HH 



HH 



(2.18) 



r,j{t+dt) = 



dts 



TO 



HH 



HH 



t + 



dt 



(2.19) 



(6) More accurate value of s{t -\- dt) and s{t + dt) is obtained using the new values of 
position, momenta and acceleration. These values are again used in step (5) to yield 
more accurate position, momenta and acceleration of the particles. These iterations 
are terminated once the values converge. 

(7) Finally, we increment the time-step and go back to step (1). 
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3. Simulation 

Molecular dynamics simulation is performed on a linear polymeric system made up of 
polymer chains. The starting point of the simulation is to generate the initial equilibrated 
melt of the polymer chains. The cell is divided into a grid of uniformly spaced points and 
the center of mass of each polymer chain is placed at these grid points. Chain growth 
takes place by generating a series of random bond vectors. A random bond vector in 
three dimensions is generated by using a random dihedral angle between to 360° while 
keeping the bond angle and bond length of the chain at their equilibrium values. This 
helps in generating an initial structure of the polymeric system very close to the one 
corresponding to the minimum energy configuration. Position of the new monomer is 
obtained by adding these bond vectors to the previously generated sites. While chain 
growth takes place, a check is performed for the occupancy of the newly generated site. 
In other words, if the new generated position for the monomer is already occupied we 
go back one step and generate a new site position. A partial overlap of the monomers 
is allowed during the chain growth which is later removed by subjecting the polymeric 
system to a soft repulsive potential to ensure that all the monomers are separated by a 
distance confirming the absence of overlap. 1 shows the initial structure of the polymer 
system generated by the above described method. 




Figure 1: Initial configuration 



Once the initial structure is ready, we perform the molecular dynamic simulation by 
subjecting the system to the forces and moments due to bonded and non-bonded poten- 
tials which describe the inter-atomic interaction. To model the van dar Walls forces, we 
use the shifted and truncated Lennard- Jones potential which acts between all pairs of 
the united atoms 



ULj{r) = 



4e 



LJ 



0, 



if r < Tc = 2.5(7 
otherwise. 



(3.1) 



Here a = 0.398 nm is the Lennard- Jones length parameter and e = 0.477 kJ/mol is the 
Lennard- Jones energy parameter that is related to the well-depth of the Lennard- Jones 
potential. Other kinds of bonded interaction used are bond stretching, bond bending and 
bond torsional potentials. 

Stretching of the bond is modeled through the finitely extended nonlinear elastic 
(FENE) potential 



UpENEir) 



1 



-KR^ In 



1 



r 

Rq 



(3.2) 



Here K = 30e/(7^ is the stiffness of the FENE potential. Rq = 1.5(7. This potential 
acts between the monomers which are directly bonded to each other in a chain. 
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i + 2 




Figure 2: Chain structure 



The bond bending potential is a function of the bond angle 9 such that 

Ubending{6) = ^KgicOsO- COS 6*0)^, (3.3) 

where Kg = 520 kJ / mol is the bond bending stiffness and the equilibrium bond bending 
angle = 110°. The bond angle 9 is defined as the angle between the two contiguous 
covalent bonds 



cos 9i 



bj 1 • 



|b,_i||bi| ' 

where b, = — rj_i is the bond vector between the atoms i — 1 and i. 
The torsional potential is a four body potential which is a harmonic function of the 
dihedral angle cj). That is 



Uu 



i{4>) = ^ A'^icos' 



(3.4) 



(=0 



where K^q = UA77 kJ/mol, K^i = -37.594 kJ/mol, K^2 = 6.493 kJ/mol and K^s = 
58.499 kJ/mol. The dihedral angle is defined as the angle between the plane formed by 
first two bond vectors and plane formed by the last two bond vectors from the set of 
three consecutive bonds 



COS( 



(b,_i X hi) ■ (b, X b,+i) 
|b,,_i X b,||bj X bi+i| 

During the simulation, the following non-dimensional parameters are used which are 
described below. 



t / € \ ■ 
time * = ^ — 

CT Vm/ 



number density p = p 



pressure p — 



temperature T — 



ksT 
e 



Table 1: Scaled parameters 

In Table 1, to is the mass of the united atom — in this case 14 amu. a and e denote 
length and energy scale parameters, respectively, and their values are the same as for the 
Lennard- Jones potential. From here onwards, we drop the bar used for scaled symbols, 
and symbols without bar should be considered as scaled parameters. In order to study 
variations and their effect on the system response, we choose a reference system with 
N = 128 polymer chains, each chain having n = 100 united atoms. The temperature of 
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this reference system is maintained at T — 4.0 in a thermal bath. The simulation cell has 
a volume consistent with a number density p = 1.0. This reference model was strained at 
a strain rate of e = 1.0 x 10^/sec. In the present simulation we use a time-step of 0.001 
for the numerical integration of equations of motion which is sufficient for the stability 
of the integration scheme and maintaining the accuracy of the same. In order to simulate 
the constant strain rate uniaxial deformation, the dimension of the simulation cell in 
the ei direction - Li - was increased at a fixed rate. We assume that the polymer is 
incompressible and apply the appropriate contraction in the transverse direction of the 
simulation cell. The polymeric system considered by us is free of cross-linkages between 
chains. The position of the atoms are also scaled in the same proportion as the dimension 
of the simulation cell. The average of the property under consideration is taken over 200 
steps of the simulation. 



4. Stress-strain and micro-structure 

First we subject the polymer system to a constant rate loading and subsequent un- 
loading at the same rate - a triangular loading - to validate our MD code with the results 
of Bergstrom & Boyce (2001). The results obtained from our MD simulation show simi- 
lar qualitative trend with those presented in Bergstrom and Boyce's study. Quantitative 
difference in the results are due to difference in the choice of model parameters, such as 
spring constants and mean bond-angles of the bond potentials, used in the respective 
simulations. 

The stress response of the polymer is correlated with its micro-structural properties 
such as bond length, bond angle, dihedral angle, mean-square bond length, mean-square 
end-to-end length, and radius of gyration distributions. The effect of the constant strain- 
rate loading on the mass ratio, mean bond angle and mean chain angle is also discussed. 



N=128, n=100, T=4.0,p=1.0, dE/dt=1.0x10^/sec N=128, n=100, T=4.0,p=1.0, dE/dt=1.0x10^/sec 
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(a) Stress-strain curve (b) Modulus vs. strain 

Figure 3: Stress response under uniaxial constant strain rate loading 

Figure 3a shows the variation of the axial stress as the system evolves in time. It is 
observed that the axial stress increases as a nonlinear function of the imposed strain 
loading. A curve fit of the stress-strain data obtained from the MD simulation is also 
shown in Figure 3a. From this curve fit we calculate the variation of the tangent modulus 
with strain. This is shown in Figure 3b. The tangent modulus initially starts with a very 
high value for small values of strain and decreases to almost zero with further increase 
in strain. In the second phase there is not much change in the modulus. However for 
strains greater than 90%, the modulus increases at a high rate. This increment is due 
to the rapid change in the mean-square bond length and mean bond angle as shown in 
Figure 4a and 4c, respectively. 

Initially the chains take the folded configuration at equilibrium. When this system 
is subjected to a constant strain rate, the linear polymer deforms due to unfolding of 
the chains as well as stretching of the bonds. Hence the tangent modulus starts with a 
high value and decreases as more and more number of chains unfold. During the second 
phase, the deformation occurs mainly due to the unfolding of the chains, and hence in 
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this phase there is very httle change in the modulus. In the third phase, chains almost 
align themselves in the stretch direction and hence any further strain directly results 
in bond stretching which is a very stiff mode of deformation. Therefore, we observe the 
increase in the modulus at high values of strain. 
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Figure 4: Variation of structural properties for uniaxial constant strain rate loading 

We now present results on the variation of the micro-structure parameters such as 
mean-square bond length, radius of gyration, mean-square end-to-end bond length, and 
mass ratio of the polymer as a function of the time-history of the constant strain rate 
loading on the polymer. 

Figure 4a shows the variation of the mean-square bond length with strain. The mean- 
square bond length gives an estimate of the contribution to the deformation coming from 
bond deformation. We observe that the mean bond length decreases initially for very 
small strains but subsequently increases uniformly with strain till approximately 60% 
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strain. After this, the bond length remains constant till 80% strain and subsequently 
increases rapidly between 80 — 100% strain. If the chains are highly oriented in the 
direction of the loading, the predominant mechanisms of deformation are bond extension 
and bond conformation Simocs et al. (2004). This is the primary reason for the increase 
in the modulus beyond a strain value of 90%. This shows that the deformation of the 
bond is strongly correlated with the stress. A sudden change in the bond length is due 
to the alignment of the chains in the loading direction. The change in the bond length 
contributed partially to the deformation, and thereby the stress, in the polymer. 

Figure 4b shows the behavior of the end-to-end length when the polymer is stretched. 
The end-to-end bond length estimates the degree of coiling of the polymer chain. If the 
chain is fully uncoiled, the end-to-end bond length will approach a value very close to 
the mean bond length times the number of bonds. In general, its value is much lower 
than this except at highly strained conditions. We observe that initially the end-to-end 
length grows linearly with strain, but in the later part of the loading, the rate of change 
of the end-to-end length grows non-linearly. This is due to the fact that in the later part 
of the loading many of the chains are open to an extent that they do not offer much 
resistance to the further opening of the chains. Hence a large and nonlinear increase in 
the mean-square end-to-end length. 

The variation of radius of gyration with strain is shown in Figure 4c. The radius of 
gyration is a measure of the coiling of the polymer chain. If the chain is highly coiled, the 
radius of gyration will take a low value, whereas if it is uncoiled, its value will be higher. 
During the initial phase of loading, the change in the value of the radius of gyration is 
less, but it slowly grows in the later part of the loading as resistance offered for uncoiling 
reduces as the chains open up. But if we keep straining the material, a stage is reached 
when most of the chains are completely uncoiled and there will not be any significant 
change in the end-to-end length and radius of gyration. 

In Figure 4d, the variation of the mass ratio to strain loading is shown. The mass 
ratio is indicative of the shape of the chains during the deformation. The mass ratio is a 
parameter that gives a rough estimate of the orientation of the chain and distribution of 
the mass in the chain. If the mass ratio is different from unity, the distribution is non- 
spherical. We find that as the system is strained, the mass ratios in both the transverse 
directions to the strain loading, namely directions 2 and 3, decrease. Moreover, in direc- 
tion 3, the distribution of the mass is very low. Hence, one can conclude that the mass 
of the polymer is distributed more or less in a plane. It means that the chains become 
more and more like a flattened cigar when it is strained. 

Now we take a look at the variation of the mean bond angle as the system is stretched. 
This is shown in Figure 4c. We find that the bond angle, during the initial part of the 
loading, increases very rapidly. Subsequent to this, the rate of change decreases, and we 
find that from e — 0.6 to e = 0.8, there is no change in the mean bond angle. Towards the 
end of the loading, that is beyond e = 0.8, there is again a sudden increase in the bond 
angle. From Equation (3.3) this results in an increase in the bond bending potential. 
Consequently, this leads to an increase in the stress towards the end of the loading, as 
observed in Figure 3a. 

Figure 4f shows the variation of the mean chain angle. The chain angle is the angle 
between the loading direction and the end-to-end vector of the chain. We observe that 
the chain angle decreases uniformly, indicative of alignment of the chains in the loading 
direction. 

In the discussion above, we presented the variation of the micro-structure parameters as 
a function of strain. We now present variation of some of these micro-structure properties 
along the chain length. In order do so, at a given segment on the chain, we take the 
average value of the structural property over all the chains. Subsequently, we take the 
time-average of this ensemble average. We repeat this for other segments on the polymer 
chain. For instance, the mean bond length along the length of the chain is obtained by 
taking the time average of the average value of the bond-length of corresponding bonds 
in all the chains. 

We observe in Figure 5a, that the bond length has fluctuations along the chain. But 
these fluctuations are uniform throughout the chain as opposed to the observation of 
Saitta & Klein (1999) wherein they report a jump in the bond length at the ends of the 
chain. This difference is due to the fact that Saitta and Klein's study considers a single 
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Figure 5: Distribution of structural properties along the length of the chain for constant 
strain rate loading 



short polymer chain loaded by an external force. As reported in their study, the mean 
value of the bond length along the length of the chain is higher than the equilibrium 
bond length value indicative of chains being in a stretched state. 

The bond angle variation along the length of the chain is shown in Figure 5b. We 
observe that the bond angles do not deviate much from the mean value and are equally 
scattered on both sides from the mean value. The variation of the dihedral angle along 
the length of the chain, for the constant strain rate loading, is shown in Figure 5c. The set 
of dihedral angles can be divided into four groups, namely {0, 113°, 247°, 360°}. This set 
of dihedral angles correspond to the minimum potential energy of the polymeric system. 
We observe that the number of bonds falling in each segment is almost the same. 

Figure 5d shows the variation of the kinetic energy during the simulation. Throughout 
this simulation, the temperature was kept constant at a non-dimensional temperature 
T = 4.0. As a result of which the kinetic energy takes a value corresponding to the 
specified temperature. In the course of the simulation, the temperature keeps fluctuating 
about this mean value with a very small standard deviation. 



5. Control parameters, stress-strain, and micro-structure 

The control variables for the polymeric system subject to dynamic strain are the pa- 
rameters that are externally set or imposed such as strain rate, density, temperature, 
and chain length. These parameters control the behavior of the polymeric system. On 
a macroscopic scale, they affect the stress-strain characteristics. Here we study the ef- 
fect of the control parameters on the stress response and correlate it with the polymer 
micro-structure response. Not only can we now explain well known mechanical behavior 
of polymers and elastomer Wincman & Rajagopal (2000), but certain anomalous behav- 
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ior can be explained too. As before, we consider a constant strain rate loading imposed 
on the polymer. 

5.1. Strain rate 

Figure 6a shows the effect of different strain rates on the tensile behavior of the polymer. 
For low values of strain, the stress is the same for all strain rates. Only for strains e > 0.2 
that we see an appreciable effect of strain rate on the stress response. Note however that 
when the strains reach a high value, the stress almost reaches a constant value. This 
is true across all the strain rates. Figure 6b is obtained by taking the derivative of the 
above curve with respect to strain. From this figure it is clear that the material is stiffer 
at higher strain rates. Now if we consider any one of these curves, there is a plateau 
region where stiffness does not vary much with strain. This is because the deformation 
in this phase is mostly due to the uncoiling of the chain. At large values of strain, the 
modulus decreases by a small amount that could be due to structural relaxation. We now 
show below how the loading rate influences the evolution of the micro-structure. 
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Figure 6: Effect of strain rate on stress response 

We first consider the influence of strain rate on the evolution of the mean-square bond 
length with strain. This is shown in 7. Mean-square bond length grows with strain for 
all strain rates and we find that at higher rates of loading the deformation in the bond 
length is more. Also note that during the initial part of loading, the bond length does 
not increase significantly. But after a strain value of e = 0.2, it increases uniformly for 
all strain rates. 
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Figure 7: Effect of loading rate on mean-square bond length 



Effect of the strain rate on the mean-square end-to-end length for constant strain rate 
loading is shown in 8. We observe that the mean-square end-to-end length increases 
in a parabolic manner. Also, as opposed to the case of mean-square bond length, we 
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find that the mean-square end-to-end length is lower at high strain rates. This indicates 
that at high strain rates internal deformation of the chains is dominant in the over-all 
deformation of the chains. 
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Figure 8: Effect of loading rate on end-to-end length 

The variation of the mean bond angle when the system is subjected to uniaxial constant 
strain rate loading is shown in 9 at various strain rates. We observe that the change in 
mean bond angle is more at strain rate e = 2.5 x 10^/sec whereas for other strain rates 
this change is almost same. Also, the variation of the mean bond angle with strain is 
very similar for all strain rates except e = 2.5 x 10^/sec. 
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Figure 9: Effect of loading rate on mean bond angle 

10 shows the effect of strain rate on chain angle for constant strain rate uniaxial 
loading. We find that, for all strain rates, the chain angle decreases at almost constant 
rate. Furthermore, we notice that alignment of the chain in the loading direction is more 
at e = 2.5 x 10^/sec where as this alignment is similar at other strain rates. 

The variation of mean-square radius of gyration for constant strain rate loading at 
various strain rates is shown in 11. We observe that the mean-square radius of gyration 
increases with strain in a nonlinear sense. There is no significant change in the variation 
of radius of gyration with strain at different strain rates. 

12 shows the variation of the mass ratios for various strain rates under constant strain 
rate loading. Mass ratios initially decrease very slowly but subsequently decrease with 
a larger rate. Beyond a strain value of e = 0.2 we observe that the rate of change of 
the mass ratios with strain is almost constant. With increasing strain in the system, we 
observe that the mass distributions decrease linearly. This is indicative of the fact that 
chains align themselves in one direction. Further, the mass ratio gs/gi is very small which 
points to the fact that in the third direction the monomer distribution is very low. Hence, 
chains become like a flattened cigar and the distribution of the monomers is primarily in 
a plane. The rate of change of mass ratio 52/51 is low at low strain rates indicating that 
alignment of the chains in the loading direction is low at low strain rates. 
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Figure 10: Effect of loading rate on mean chain angle 
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Figure 11: Effect of loading rate on radius of gyration 
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Figure 12: Effect of strain rate on mass ratios 

We find that at a high value of strain rate the change in the mean-square bond length 
and mean bond angle is more, where-as the change in the mean-square end-to-end length 
is small and mean-square radius of gyration remains same at all strain rates. This indi- 
cates that internal deformation in the chain dominates over the uncoiling at higher strain 
rates. 

5.2. Density 

The effect of density on the stress-strain behavior of the polymer for a constant strain 
rate is shown in Figure 13a. Higher the density, greater is the stress in the polymer. Also, 
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we find a very huge jump in the stress from p = 0.8 to p = 1.5. In fact, jump in the stress 
level from p = 1.2 to p = 1.5 itself is very significant. Similar behavior was also observed 
by Bower & Weiner (2006). Bower & Wcincr (2004) report glass transition takes place 
at p = 1.2 and T = 2.0 in which they mention about presence of an additional energy 
component in the difference stress. This component increases with increment in density 
and decrement in temperature. This could be one of the possible reasons for a huge jump 
in stress from p = 1.2 to p = 1.5. 

Variation of the modulus is shown in Figure 13b. One observes three phases of defor- 
mation in the material. During the initial loading phase the polymer with higher density 
is more stiffcr. In the intermediate range of strain values, for the polymer with density 
p = 1.5, there is an increase in stress and even the modulus increases too. This is due 
to the fact that for a dense polymer it takes more time to relax. This is not true in 
the modulus versus strain curves for lower values of strain. This intermediate or plateau 
region is followed by a third region where the modulus again increases. This is due to 
the chains being sufficiently uncoiled and the imposed strain resulting in the deforma- 
tion of bonds that are stiff. Therefore, the stress and the modulus in the polymer starts 
increasing during this phase. 
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Figure 13: Effect of density on stress response 
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Figure 14: Effect of density on mean-square bond length 



We now study the effect of the density on the mean-square bond length as the polymer 
is strained. The effect is shown in 14. There is a critical value of density, p — 1.5, above 
which mean-square bond length increases with strain. At values of densities below this, 
the mean-square bond length shows a rather slow increase with strain. This implies that 
deformation of the bonds at higher densities is more. Also, at high density values, the 
chain is tightly packed and sliding and uncoiling are not easy because of the constraints 
from the neighboring molecules. As a result, bond deformation is the only mechanism to 
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provide the appropriate strain in the material. When there is no strain in the system, we 
observe that the mean-square bond length takes a lower value at higher densities to give 
a compact configuration. 
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Figure 15: Effect of density on end-to-end length 



The effect of density on the mean-square end-to-end length was also studied for con- 
stant strain rate loading. The results are shown in 15. At a lower value of density, the 
mean-square end-to-end length is larger than that at a higher value of density. This is 
since at lower densities the chains get enough space to uncoil themselves when subjected 
to strain. This is not possible at higher densities. Also, we observe there is not much 
difference in the mean-square end-to-end length beyond p = 1.0. The only significant 
difference is observed for a density value of p = 0.8 in which case the chains have enough 
space in between them to open up. 
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Figure 16: Effect of density on radius of gyration 



16 shows the effect of density on the radius of gyration under the constant strain 
rate loading. At lower densities, the radius of gyration is larger than that at higher 
density for the same reason as with the mean-square end-to-end length discussed above. 
But, this effect is more dominant for density p — 0.8. For other values of densities, this 
effect is not very significant. The radius of gyration initially increases slowly with strain 
but subsequently increases at a faster rate because in this phase the chain uncoiling is 
predominant. 

The variation of the mass ratios under constant strain rate loading at different densities 
is shown in 17. Higher the density, the initial configurations of the molecules are more flat 
to provide the compactness. This remains true approximately at all the values of strain 
during stretching. As the system is stretched, initially the change in the mass ratios is 
very small but subsequently decreases with constant rate. At a density of p = 1.5, this 
transition from a small change in the mass ratios with strain to a change at a constant 
rate is delayed and we see the cross over of this curve with the p = 1.2 curve. Mass ratios 
decrease with strain due to alignment of the chains in one direction. 
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Figure 17: Effect of density on mass ratios 
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Figure 18: Effect of density on mean bond angle 



The variation of the mean bond angle at different densities for uniaxial constant strain 
loading is shown in 18. We find that before the loading is applied, the mean bond angle 
is lower at higher densities to give a compact configuration. As the system is stretched, 
mean bond angle increases at all the densities. But, the rate at which the mean bond 
angle increases is greater at higher values of densities. At low densities, the change in the 
mean bond angle is very small. Note that irrespective of the density of the polymer, in 
all the cases, the mean bond angle converges to the same value at high values of strain. 
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Figure 19: Effect of density on mean chain angle 
The effect of density on mean chain angle for uniaxial constant strain rate loading is 
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shown in 19. There is no significant effect of density on the behavior of mean chain angle 
and it decreases with loading in similar manner at all densities. 

5.3. Temperature 

The stress response of the system at different temperatures is shown in Figure 20a. One 
observes that lower the temperature, higher is the stress in the polymer. This was also 
observed by Clhui & Boyce (1999) who performed a Monte-Carlo simulation of a polymer 
in a compression test. At lower values of temperature, the kinetic energy of the united 
atoms is very small and hence it takes more time for uncoiling and relaxation. Whereas 
at higher temperatures, the monomers have high kinetic energy and hence relax faster. 
Therefore, we find the stress developed in the polymer at higher temperatures is less as 
it attains a more relaxed structure very soon. Also, at higher temperatures, distribution 
of the monomers in the chains change from spherical to a less spherical shape slowly 
under the stretch. In contrast, at low temperatures, this change is slow initially but 
subsequently becomes faster. This is also corroborated by the mass ratio variation with 
strain at different temperatures shown in 24. This indicates that at high temperatures 
anisotropy in the stress is low resulting in lower stress. The modulus curve for different 
temperatures is shown in Figure 20b. This figure clearly shows that at low temperature 
the polymer exhibits a higher modulus at higher strains but at low strains there is very 
small difference in modulus. 
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Figure 21: Effect of temperature on mean-square bond length 



We begin our investigation of the influence of temperature on the micro-structure 
parameters that determine stress response with the effect of temperature on the mean- 
square bond length. The result is shown in 21. We find that the mean-square bond 
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length is larger at higher temperatures. Also, we note that as the system is loaded, the 
mean-square bond length at different temperature values increases with the same rate. 
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Figure 22: Effect of temperature on end-to-end length 



The effect of temperature on the mean-square end-to-end length under constant strain 
rate loading is shown in 22. At higher temperatures the end-to-end length is larger 
because the kinetic energy dominates over the potential energy and the system relaxes 
faster. 

However, note that though the mean-square bond length and the mean-square end-to- 
end length variations with strain show significant change at different temperatures, the 
axial stress versus strain variation at different temperatures, shown in Figure 20a, does 
not show any significant change. This is primarily because though the temperature does 
affect the bond length as well as the end-to-end length significantly, the orientation of the 
bonds and chains being random, stresses generated are isotropic and therefore increase 
the pressure in the polymer leaving the anisotropic axial stress largely unaffected. 
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Figure 23: Effect of temperature on radius of gyration 



The effect of temperature on mean-square radius of gyration in constant strain rate 
loading is shown in 23. We observe that there is no significant effect of temperature on 
mean-square radius of gyration and its variation remains similar. This is also reflected 
in the stress-strain behavior where we do not find any significant variation at different 
temperatures. 

The variation of mass ratio under uniaxial constant strain loading at various temper- 
atures is shown in 24. There is no significant difference in the variation in the mass ratio 
with temperature except that it decreases under tensile strain at all temperatures. We 
observe that at high temperatures, chains are more flat initially but as the system gets 
stretched, initially change in mass ratio (72/31 is very slow at higher temperatures and 
it continues till higher values of strain. Subsequently, mass ratio at all the temperatures 
decreases at constant rate. As a result we find the cross-over of the variation with 
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Figure 24: Effect of temperature on mass ratios 

strain at different temperatures, g-ij g\ is lower at higher temperatures, the only exception 
being at T = 5.0 where this curve falls between T = 3.0 and T = 4.0. 
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Figure 25: Effect of temperature on mean bond angle 

The effect of temperature on mean bond angle is shown in 25 for uniaxial constant 
strain rate loading. We observe that the mean bond angle increases with the loading 
in all the cases, but, there is no significant difference in mean bond angle at various 
temperatures. 
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Figure 26: Effect of temperature on mean chain angle 



26 compares the variation in mean chain angle at various temperatures when the 
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system is subjected to uniaxially constant strain rate load. We find that the mean chain 
angle again decreases with strain at all temperatures with more alignment observed at 
higher temperature. But, this change is again not very significant leading us to state that 
temperature does not affect the behavior of the mean chain angle. 

5.4. Chain length 

We consider different chain lengths of the polymer and subject the system to constant 
strain rate. The chain lengths are varied by varying the number of monomers in a polymer 
chain. The density is kept constant by controlling the simulation cell dimensions. The 
stress response is shown in 27. Short length polymers, n = 10, 20, and below, show 
distinctly different stress response. Longer polymers develop more stress and we observe 
that their modulus or stiffness is greater. A study by Krcmer & Grest (1990) shows 
that below a chain length = 35 there is no entanglement in the polymer chain. This 
length is called the entanglement length. Our observation of increased stress in long chain 
polymers, that is polymers whose chain length in greater than n — 50, can be attributed 
to entanglement effects. We also observe that beyond a particular value of chain length, 
n = 100, there is no significant difference in the stress-strain behavior. 
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Figure 27: Effect of chain length on stress response 
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Figure 28: Effect of chain length on mean-square bond length 



The effect of chain length on the mean-square bond length is illustrated in 28. When 
the chains are of smaller length, there is no significant change in the mean-square bond 
length and the deformation is mostly from the sliding of the chains on one another. 
Short chain polymers can slide over each other very freely as possibility of entanglement 
is rare. Also, since the chains are shorter, there is not much scope for uncoiling of chains. 
Further, we observe that at high strain values, in the range 0.9 < e < 1.0, the mean- 
square bond length decreases for shorter polymers. In fact, this is the main reason for 
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the drastic change in their stress-strain behavior as observed in 27. In the case of long 
chain polymers, the mean-square bond length grows with strain since chains cannot slide 
so freely. In such a case, the total elastic deformation has contributions from the partial 
sliding of the chains, uncoiling, and the deformation due to stretching of the bonds. 
Also, beyond n = 100 there is no significant change in the variation of the mean-square 
bond length. This is one of the reasons for similar stress-strain behavior of longer chain 
polymers as seen in 27. 

29 compares the effect of chain length on the mean-square end-to-end length under 
constant strain rate loading. The mean-square end-to-end length shows very little varia- 
tion with strain for chain lengths n = 50 and below. This is since for small chain lengths 
there is not much scope of uncoiling. For longer chain lengths, uncoiling is possible and 
in fact at higher values of strain there is a nonlinear increase in the end-to-end length. 
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Figure 29: Effect of chain length on end-to-end length 
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Figure 30: Effect of chain length on radius of gyration 

The effect of chain length on the radius of gyration is presented in 30 for constant strain 
rate loading. For shorter chains, as expected, there is no significant change in the radius 
of gyration. But, longer chains show a larger increase of radius of gyration with strain 
and also a increasing rate of change with strain. In short chain polymers, the overall 
deformation is due to the alignment of the chains in the loading direction. Further, there 
is slipping over one another rather than actual deformation of the chains. 

The effect of chain length on the mass ratios in constant strain rate loading condition 
is shown in 31. We observe that for very short chains, say chain lengths with n = 10 or 
less number of monomers, the mass ratio is very small in the beginning of the loading. 
This implies that very short chains are flattened from the beginning. For larger chain 
polymers the mass ratios take higher value in the beginning of the loading. But, beyond 
n = 100 we observe that mass ratio again take lower values as the chain length increases. 
This indicates that very longer chains are flat because of the constraints coming due 
to presence of neighboring molecules. Now, as the system is stretched, we find that in 
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Figure 31: Effect of chain length on mass ratios 
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very short chain polymers, change in mass ratio is small as compared to the that in long 
chain polymers. As a result, in a long chain polymer, the mass is distributed more in 
the loading direction at high values of strains. For n = 500, we observe that the chains 
almost take a one dimensional structure. 
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Figure 32: Effect of chain length on mean bond angle 



The effect of chain length on bond angle is shown in 32 for constant strain rate loading. 
For very short chain polymers (n = 10) we find that there is a lot of fluctuations in the 
bond angle but the overall change in it is very small as compared to longer polymers. Even 
at n = 20 we observe fluctuations at high values of strain. The variation is similar for all 
chain lengths and it increases as the system is pulled. Small chain length polymers can 
easily slide over each other and hence the deformation is mostly due to slipping rather 
than actual internal deformation, a fact noticed while discussing 28, 29 and 30. Long 
polymers cannot slip over each other so easily and encounter entanglement constraints. 
Hence, internal deformation of the chains is dominant in this case and we see uniform 
increment in the mean bond angle with strain. 

33 shows the effect of the polymer size on the mean chain angle. We find that at the 
initiation of loading itself, the mean chain angles are lower for longer polymers. As the 
system is pulled, the alignment of the chains in the loading direction is observed for all 
the polymer chain lengths considered. But note thought that the alignment is typically 
more for longer polymer chains. This is also supported by the initial alignment of the 
chains in case of long polymers. Except for n = 200 we find that as the system size 
increases, the alignment of the chains throughout the loading is more for long polymers. 
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Figure 33: Effect of chain length on mean chain angle 



6. Conclusions 

We have presented a study of the stress response of a linear polymer subjected to 
uniaxial constant strain rate loading. The tool used for the study was molecular dynamics 
simulation. The distribution of micro-structure parameters along the chain length as well 
as the time-evolution of the micro-structure parameters is computed. These are then used 
to understand the effect of the micro-structure parameters on the stress response. We 
observe that the bond length has fluctuations along the chain. But these fluctuations are 
uniform throughout the chain as opposed to the fact that a jump occurs in the bond 
length at the ends of the chain as in a single short polymer chain loaded by an external 
force. While investigating the effect of externally defined parameters such as temperature, 
density, chain length, and strain rate, we observe that density has a very strong effect on 
the stress-strain behavior of the polymer and there is a large jump in the stress levels as 
the density is increased. At higher densities the deformation is mostly dominated by bond 
deformation, and change in overall size and shape of the polymer chain is less significant. 
Characteristics of the polymer are very different for different chain lengths. Short chain 
polymers more or less behave like rigid molecules. There is no significant change in their 
internal structure when loaded. On the other hand, long chain polymers are very sensitive 
to the loading and the external environment. In the case of short chain polymers, those 
below the entanglement length, since there is no significant change the internal structure, 
the deformation is mostly due to slipping of the chains over one another. In long chain 
polymer this is not the case as they get entangled. Deformation in long chain polymers 
is due to deformation of the internal structure as well as their uncoiling. Temperature 
does not have a very significant effect on the stress strain behavior, but, it has significant 
effect on the mean-square bond length of the polymer and we find that relaxation in 
the polymer is faster at higher temperatures. Further, the end-to-end length, radius of 
gyration, and other parameters that characterize the micro-structure are not affected 
much by temperature. But, we find that chains loose their coiled structure very slowly at 
high temperature. Rate of the loading mostly affects the micro-structure very mildly and 
thereby stress-strain response. We observe that parameters associated with the covalent 
bond increase with the rate of loading where as parameters representing the overall size 
of the polymer chains decrease. This indicates that internal deformation of the chain 
dominates at higher rates of loading. 

The united atom model used for the molecular dynamics simulation is a powerful tool 
to study the statistics of polymer dynamics and correlate it to mechanical behavior. 
This can then be useful to design polymers and elastomers with desired mechanical 
characteristics. 
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